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ABSTRACT 

In this paper we study the wind accretion onto a rotating black hole in the close binary system 
harboring a young massive star. It is shown that the angular momentum of the accreted stellar 
wind material is not sufficient for the formation of an accretion disk. On the other hand, in the 
considered conditions the Blanford-Znajek mechanism can be activated, thus powerful jets 
can be launched in the direction of the rotation axis of the black hole. Importantly, no obser- 
vational signatures of accretion, as typically seen from the thermal X-ray emission from the 
accretion disks, are expected in the suggested scenario. Here, properties of the generated jet 
are studied numerically in the framework of a 2D general relativity magnetohydrodynamical 
approach. Due to the accumulation of the magnetic flux at the black hole horizon, the jet power 
is expected to be modulated on a sub- second time- scale. Although the intervals between jet 
active phases depend on the magnetic flux escape from the black hole horizon (which can be 
modeled self-consistently only using a 3D code), a general estimate of the averaged jet power 
is obtained. It is expected that for the black hole rotation, expected in stellar binary system 
(the dimensionless rotation parameter a = 0.5), approximately 10% of the accreted rest energy 
can be channeled into the jets. In the specific case of the gamma-ray binary system LS 5039, 
the obtained jet luminosity can be responsible for the observed GeV radiation if one invokes 
Doppler boosting, which can enhanced the apparent flux from the system. 
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1 INTRODUCTION 

Several binary systems, consisting of a massive star and a com- 
pact object (CO) -neutron star (NS) or black hole (BH)- were 
detected in the very high energy (VHE) regi me ([Aharonian et al. 



2005b a 2006 20091 [Albert et al.|2006|[2057l|2009[|Acciari et al. 



2008t|2009t |2011 ), and their non-thermal emission is believed to 
be linked to the processes related to the CO (see e.g. |Bosch-Ramon[ 
|& Khan gulyan 2009j for a review). Moreover, given large val- 
ues of gamma-gamma opacity expected in the close vicinity of 
CO, the VHE gamma-ray production most likely should be related 
in the outflow originated in the CO. In particular, in the system 
PSR B1259-63/LS2883 a pulsar is orbiting a massive 0-type star 
(see [Neguemela et al.|201 Ij and reference therein), and the detected 
X-ray and TeV radiation is attributed to synchrotron and inverse 
Compton (IC) emission produced by electrons accelerated at the 
pulsar wind termination shock (Kirk et al.|2000[ [Khangulyan et al.[ 
|2007[ [Uchiyama et al.^2009 1, w hile emission detected in the GeV 
energy range (Tam et al.|20lT|[Abdo et al.|201l] ) can be explained 



by the bulk Comptonizaiton of the pulsar wind I Khangulyan et al.[ 
|2011| l. In the case of Cygnus X-1, the evidence of detection in TeV 
found MAGIC (Albert et al. 2007) is most likely related to the non- 
thermal activity occurring in the jet launched by a black hole (^Esinj 
|etal.|1998[ ) accreting matter from a 09.7Iab star ( [Zi61kowski|2005[ ). 

The classification of compact objects (CO) in these two cases 



is rather clear since (i) pulsed radio emission coherent with the or- 
bital motion is detected from PSR B1259-63/LS2883 |Johnston| 
,et al. 1996 ); and (ii) in the case of Cygnus X-1 the mass of the CO 
exceeds the Oppenheimer-Volkoff upper limit for the mass of NS 
( [Casares|2007j), and the accre tion disk is clearly seen in the X-ray 
energy band ( Esin et al.|1998| ). 

The physical processes responsible for the VHE radiation de- 
tected from two other systems, LS 5039 and LS I +63-303, is still 
debated, since currently there is no observational evidences which 
allows robust conclusions regarding the nature of the CO. In par- 
ticular, despite intensive searches no pulsed emission has been de- 
tected from these systems. Although, strong free-free absorption, 
expected in these relatively compact binaries, could hide the radio 
pulsations ( |Dubus|2006|), the stri ct upper limits obtained in the X- 
ray energy band ( Rea et al.|2010[[20TT| l pose additional constraints 
for the realization of the binary pulsar scenario in these systems. 
Another, independent on the orientation of the "pulsar beam", an- 
other test for the presence of a pulsar in the binary system is re- 
lated to the impact of the pulsar wind on the stellar environment. 
Namely, the collision of the stellar and relativistic pulsar winds 
should give rise not only to the non-thermal emission related to 
the pulsar wind, but the stellar wind would be also shocked pro- 
ducing thermal X-rays. This approach was used to constrain the 
spin-down (SD) luminosity of the possible pulsar in the case of 
LS 5039 ( [Zabalza et al.|201l] ). Importantly, the derived upper lim- 
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its appeared to be very close to the non-thermal luminosity detected 
with Fermi Large Area Telescope (Fermi/hAI) ( |Abdo et al.|2009a| ). 
Thus, one should expect either a perfect conversion of the pulsar 
energy to GeV gamma-ray, or to invoke Doppler boosting of the 
emission. Although, the relativistic motion of the shocked pulsar 
wind is a rather natural effect in binary pulsar systems ( Bogov-| 



|alov et al.|2008l|Khangulyan et al. 12008) |Dubus et al.|20iot , fast 
mixing of the stellar and pulsar winds may destroy the relativis- 
tic regime of the outflow (Bosch-Ramo n & BarkoV|2011j ), and the 
efficient Doppler boosting of the emission is possible only on the 
scales limited by the binary system size. 

Orbital phase dependent extended radio structures, predicted 
for the binary pulsar systems Dubus (2006), were found in the 
gamma-ray bin ary systems PSR B1259-63/LS2883, LS 5039 and 
LS I +63 -303 ( IDhawan et all2006l|Rib6 et al.|2008[|Mold6n et al.| 
|2011a|b| ). This may be considered as a strong evidence for the pres- 
ence of pulsars in these systems, however (i) the extended radio 
emission can be naturally produced in gamma-ray binary systems 
via the synchrotron emission of secondary pairs ( )Bosch-Ramon| 
|et al.||2008| [Bosch-Ramon & Khangulyan||201ll ); (ii) to interpret 
the resolved structures requires a knowledge of the flow dynam- 
ics that presently lacks. This uncertainty can be to some extend 
explained by the excessive simplification adopted in the theoreti- 
cal calculations, in particular the interaction of the winds in binary 
pulsar can be rather complicated ( [Bogovalov et al.|[20 08 , Bosch- 
Ramon & Barkov|2()Tl ) and diff'er significantly |Bogovalov et al.| 



from the case realized at the interaction of a pulsar with the 



interstellar medium ( [Kennel & Coroniti|1984l|Dubus|200'6t . Thus^ 
the interpretation of these observations in the context of the pulsar 
scenario remains inconclusive. 

An alternative scenario implies the production of VHE 
gamma-rays in miroquasar jets launched by an accreting black 
hole. Although both systems, LS 5039 and LS I +63-303, do not 
show thermal emission produced by the accreation disk, as it is 
expected in the classical microquasar (piQ) paradigm, the spectro- 
scopic observations of LS 5039 favor a mass of the CO (Casares 
|et al. 12005) [Sarty et al.|2011| exceeding the Oppenheimer- Volkoff" 
upper limit for the mass of NS. If a BH is indeed located in these 
systems, the absence of the thermal emission of the accretion disk 
might be interpreted as a signature of the radiatively non-efficient 
disk accretion ( [Narayan & Yi|1994||BTaridford & Begelman|1999| ), 
although the applicability of these scenarios is debated for sub- 
Eddington stellar systems ( [Artemova et al.|2006 |. Therefore, a de- 
tailed study of the properties of the accretion flow in binary system 
is required for meaningful conclusions regarding the possibility to 
power these sources, LS 5039 and LS I +63-303, by accretion. We 
show below that in these systems of massive stars, which do not fill 
the Roshe lobe and the CO is not very massive, the accretion onto 
the BH should occur from the stellar wind. 

For the sake of specificity, we adopt in this work the system 
parameters similar to those of LS 5039. In this system the CO 
is orbiting an 06.5V star in a relatively non-eccentric orbit (see 
the system parameters in Tab. 1). This system was detected as a 
bright non-thermal source in a very broad energy range, from radio 
to VHE (see e.g. |Bosch-Ramon & Khangulyan[[2009[ and refer- 
ences therein). The spectral energy distribution (SED) of the radi- 
ation is dominated by GeV gamma-rays. Namely, a flux of about 
fewxlO^^ ergcm~^s~^ was detected with Fermi/LAI (Abdo et al. 



|2009b| ). The fluxes in X-ray and TeV ener gy bands are compara 
ble and at the level of ~ 10^"^ ergcm~^s~^ ( [Aharonian et al.|2006 



[nian et al.[20'06)[Kishishita et al.|2009[[Abdo et al.|2"009b) ). We note 
that the emission detected from LS I +63-303 is known to diflfer 
significantly from orbit to orbit ( |Li et al.|2011)[Acciari et al.|2()TT| l. 

According to our analytical estimates, the accretion rate can be 
relatively high, in agreement with the results obtained by [Owocki| 
|et al.[ ( [20TTi| ). Importantly, due to the lack of angular momentum, 
the accretion proceeds without the formation of an accretion disk, 
i.e. one should not expect detectable thermal X-rays from the in- 
falling matter. Despite the absence of the accretion disk, our study 
shows that a powerful jet can be launched in this accretion regime. 
The crucial conditions for the formation of a powerful jet are (i) 
magnetization of the in-falling flow; and (ii) high enough accre- 
tion rate. If these conditions are fulfilled, a fast rotating BH can 
power and launch a magnetically driven jet ( [Komissarov & Barkov| 
[2009 1 ). Accretion of the matter leads to the accumulation of the 
magnetic flux near the BH. When the B -field reaches the critical 
value, the magnetic pressure pushes out the in-falling matter paus- 
ing the accretion and dis-power the jet. This phenomenon is similar 



& Ruzmaikin|1974[|1976)[Kaisig et al.|1992 


Narayan et al.|2003| 


Igumenshchev[[2008| Romanova et al.[f2009 


Tchekhovskoy et al.| 



2011 ). The periodic stops of the central engine leads to a forma- 
tion of a quasi-periodic structure of the jet on a scale significantly 
smaller (by three orders of magnitude) than the system semi-major 
axis. 

The paper is organized as following: in Sect[2] we present 
the general properties of the accretion flow in binary systems; in 
Sect. |2.2[ we develop a 2D general relativity magnetohydrodynam- 
ical (GRMHD) numerical model for the accretion flow, and study 
difl'erent accretion regimes; in Sect. |3] we discuss the implications 
of the obtained solutions; and in Sect. |4] we summarized the ob- 
tained results. 



2 MODELING ACCRETION 

2.1 General Properties of the Accretion Flow 

In a binary system harboring the BH, and in which the optical com- 
ponent does not fill the Roshe lobe, the BH will accrete from the 
stellar wind directly. The accretion rate can be estimated from the 
Bondi-Hoyle solution, which describes the accretion onto a cen- 
ter of mass M moving with a constant speed Voo through a uniform 
medium of density poo- The Bondi-Hoyle mass accretion rate is well 
approximated by the equation 



Mbh = nR\p ooVo. 
where 

2GM 



(1) 



(2) 



is the accretion radius. In the Bondi-Hoyle solution the mean angu- 
lar momentum of accreted matter is zero. 

The wind accretion can be only marginally described by the 
Bondi-Hoyle solution due to the density and velocity gradients 
across the BH trajectory, although provided that the gradients are 
small on the scale of the a ccretion radius, the abo ve expression 
for M is still quite accurate ( Ishii et al.|l993 



Ruff^ert 1997 1999). 



[Takahashi et al.|2009] ). X-ray, GeV and TeV emission components 
were shown to be periodic and stable on year time- scales { |Aharo-[ 



However, the accreted matter can obtain a non-zero mean angular 
momentum flllarionov & Sunyaev'1 975| [Ishii et al.|1993[ [Davies[ 
& Pringle 1980, Rufl^ert 1997 1999|71n analogy to the work by 
Barkov & Komissarov| ( |201 1) given the unsettled nature of this is- 
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Table 1. The parameters of the system LS5039 



Description 


Designation 


Value 


Mass of star 


Ms 


26Mo 


Radius of star 




93Rq 


Temperature of the star 


Ts 


39, 000 K 


Stellar Wind termination velocity 


Voo 


2, 400 km/s 


Stellar Wind loss rate 


Ms 


4 X lO-^Moyr"^ 


Orbital period 


Ps 


3.9 day 


Eccentricity of the orbit 


e 


0.24 


The mass of the BH 


Mbh 


3Mo 


Semimajor axis 




3.5/?, 



sue, we will assume that the mean specific angular momentum of 
accreted matter inside the accretion cylinder is 

<A> = ^nRl, (3) 

where Q is the angular velocity of the orbital motion; and 77 is a free 
parameter (\t]\ < 1), which reflects our current ignorance. 

In the specific case of LS 5039, the above estimates allow us 
to obtain analytically approximate values for the accretion rate and 
angular momentum of the flow. The physical properties of the sys- 
tem are summarized in Table[T](see |Casares et al.|2005[|Sarty et al.| 
|2011 . and the references therein). 

The wind velocity profile can be approximated as ( |Kudritzki| 
|& Puls|2000| ) 

v^{R)^v^{l-^y (4) 

Thus, the typical wind velocity at the orbital separation distance, ^o, 
can be estimates as UwC^^) ~ 1.7 x 10^ cm s~^ (see Table. 1, for the 
used parameter values). This value exceeds significantly the orbital 
velocity of the BH, ~ 4 x 10^ cm s~\ thus in the calculation one 
can safely neglect the BH velocity. Substituting the wind speed to 
Eq.([TJ one obtains the value for the BH accretion rate: 

Mbh = 1.5 X IQ-^^Moyr"^ . (5) 

We note that this value depends on the star separation distance and 
can vary up to a factor of 3-5 along the orbit ( Owocki et al. 2011 ) . 
Note, that the accretion rate in the system LS 5039 is significantly 
smaller than in the case of Cygnus X- 1 due to large diff'erence in 
the masses of CO and in the speed of the stellar wind. 

Follow Eq.|2]), which gets Ra ^ 2.7 X 10^^ cm, after it substi- 
tution in Eq.([3]) one can obtain the value for the angular momen- 
tum of the accreted matter Ja = 3.3 x 10^^ cm^ s~^ This value 
appears to be small as compared to the specific angular momen- 
tum related to the last marginally bound orbit, jmb ~ 3.4r^c = 
3AGMbh/c - 4.5 X 10^^ cm^ s~^ (here we assume the BH dimen- 
sionless rotation parameter to be a = 0.5) ( [Bardeen et al.||1972| ). 
This relation, <^ jmb, implies that the accretion disk cannot 
be formed in such an accretion regime, and one should expect a 
nearly spherically symmetric in-flow of the stellar material. This 
allows to perform the calculations of the accretion under 2D ap- 
proximation. Namely, we study the case of the spherical accre- 
tion onto a rotating black hole (the setup is similar to ones consid- 
ered by Komissarov & Barkov|2009[[Barkov & K omissarov |2010| ). 
The initial magnetic field is assumed to be homogeneous on the 
scale of the computational domain Rcom - 2 x 10^ cm, and par- 
allel to the rotation axis of BH. The accretion rate was assumed 
to be steady on the computational time-scale, and two values were 
adopted Mbh = (1-5 and 8) x IO-^^Mq yv'^ . 



2.2 Numerical results 

The main details of our numerical method and various test simu- 
lations are described in the literature (see for details |Komissarov| 
[T999I [20041 |2006l [Komissarov & Barkov|2007 V Here we outline 
just the key elements of the approach. The calculations are per- 
formed in Kerr-Schild spacetime coordinates. The computational 
grid was selected to be uniform in polar angle, 6, where it has 
320 cells and logarithmic in spherical radius, r, where it has 773 
cells. The inner boundary is located just inside the event horizon 
and adopts the free-flow boundary conditions. The outer bound- 
ary is located at r = 2.2 x 10^ cm and at this boundary the flow 
is prescribed according to the free-fall Bondi model with constant 
accretion rate. 

In our study we check the feasibility the Blanford-Znajek 
(BZ) driven jet form ation (Ruffini & Wilson^ 1975) [Lovelace] 1976[ 
jBlandford & Znajek[[1977| ) in the specific conditions of wind ac- 
cretion in the gamma-ray binary system LS 5039. Simple relations 
obtained by [Barkov & Komissarov ( 2008bj relate the BZ energy 
realize to the magnetic flux 

Lbz = 1.4x10'V(«)^2o(^) ergs-', (6) 

where f{a) - (\ + Vl - a^) is a dimensionless function ac- 
counting for the BH rotation; and ^20 = ^/lO^^ Gcm^ is the 
magnetic flux. According to [Barkov & Komissarov| ( [20T0| ), the ex- 
pected value of the dimensionless rotation parameter of the BH in 
a stellar binary system should be relatively small, namely a - 0.5 
(/(a) = 0.072). 

Regarding the strength of the initial magnetic field, to reduce 
the calculation time we have assumed a fixed value of 10^ G, which 
exceeds significantly the value of a few G expected in stellar binary 
systems. We note that the actual value of the magnetic field strength 
aff'ects only the duration of the magnetic flux accumulation stage, 
but does not influence the energy release at the saturation phase (see 
details below). Thus, the assumption of a strong magnetic field does 
not influence our results while allowing to reduce significantly the 
computation time. 

At the initial stage of the simulation spherically symmetric ac- 
cretion occurs and the magnetic flux accumulates on the horizon of 
the BH. At the moment when the magnetic flux exceeds the criti- 
cal value (see [Barkov & Komissarov|2008a|[Komissarov & Barkovj 
|2009| for more details), the BZ mechanism gets activated and a 
magnetically driven jet is launched. The accumulation phase can 
be clearly seen in Fig.[T] as the lower trajectory in the magnetic flux 
vs square root of the jet power phase space. Namely, during this 
period the system is characterized by an increasing magnetic flux 
without any ejection of energy. This is caused a too low value of 
the magnetic torque to activate the BZ mechanism. 

The flux accumulation stage (marked with number 1 in Fig[TJ 
continues until the magnetic field flux reaches the critical value of 
^ ^ lO^^G cm^ for an accretion rate of Mbh = 8 X lO'^^Moyr-^ 
(see transition in Fig[T] left panel). This value of magnetic flux can 
be used to check the feasibility of the jet activation in binary sys- 
tems with relatively weak magnetic field (we remind that in the 
calculations we have assume a large value of B -field for the sake of 
reducing of the computational expenses). The required value can 
be reached even in the case of accretion of an uniform magnetic 
field of strength ~ 0.1 M^f^ G from the radius /?a - 3 x 10^^ cm. 
We note that stellar winds can carry significantly stronger magnetic 
fields (see eq. [Usov & Melrose] 1992 1. 

Once the jet is launched, the dense falling plasma is evacu- 
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ated from the BH magnetosphere, and the BZ mechanism starts 
to operate in the conventional way. Thus, the jet power should be 
described by Eq. ([6]). Indeed, in Fig.[T]one can see that after the ac- 
cumulation/cleaning (marked with numbers 1 and 2 in the figure) 
stages the system remains on a straight line Lg^ ^ ^ (marked with 
number 3 in the figure). 

The jets are ejected in the directions of the BH rotation axis 
and blow up bubbles, as it is shown in Fig|2](right top panel). In the 
zoomed image of the central part, shown in Fig |2] (left top panel), 
one can see that the matter in-falling occurs through a narrow re- 
gion in the equatorial plane. However we note, that although this 
can mimic a thin accretion disk ( Bisnovatyi-Kogan & Ruzmaikin| 
|1974| |1976| ), the matter velocity has no rotational component in 
this region. This can be clearly seen from the left bottom panel in 
Fig [2] where the ratio of toroidal and poloidal components of the 
B -field are shown. Any toroidal velocity component would lead to 
an enhancement of the corresponding (i.e. toroidal) magnetic field 
component. It is noteworthy that this configuration with radial ac- 
cretion along the equatorial plane does not significantly diff'er from 
the standard disk accretion (which was adopted in the conventional 
BZ mechanism) in the regions close to the BH horizon since, even 
in the case of disk accretion, a free falling region is formed inside 
the marginally stable orbit. 

Another important aspect of the considered accretion regime 
is shown in Fig|2] bottom right panel, where the ratio of gas and 
magnetic pressures is shown. It can be seen that in the equatorial 
region the gas pressure dominates, but the jet itself is strongly dom- 
inated by the magnetic pressure. This configuration fits perfectly 
the paradigm of the magnetically driven jets, where the gas pres- 
sure in the jet is expected to be negligible, while in the equatorial 
plane the presence of matter is required to keep a strong magnetic 
field near the horizon. 

The time dependence of the jet power is shown in Fig [3] It 
can be seen that after the magnetic flux accumulation stage, the jet 
power starts to increase quickly. This growing phase is stopped by 
a sudden drop of the jet power. This is caused by gradual increase 
of the magnetic field flux at the BH horizon during accretion. When 
the magnetic pressure overtakes the in-falling matter ram pressure, 
a fast change of the structure of the B -field occurs, leading to the 
formation of a magnetically arrested torus, as shown in Fig[4j A 
similar situation is discussed by Narayan et al., ( ,2003] ); |Igumen-| 
[shchevlpOOSI for the case of disk accretion. 

Since in the case of a magnetically arrested torus no accre- 
tion occurs, a fast increase of the gas pressure is unavoidable. Once 
the gas pressure reaches a value exceeding the B -field pressure, the 
accretion, as well as the jet power supply, are resumed. We note 
however that in the case of a 3D setup a natural escape of the mag- 
netic flux can be realized through the generation of magnetic tubes. 
In principle, this could lead to the formation of a nearly steady flow, 
although a quasi periodic pattern may still remain in the jet power 
time dependence. A detailed study of this issue is beyond the scope 
of this paper and will be discussed elsewhere, below we present 
a qualitative discussion of this efl'ect, which despite of qualitative 
still allows rather fundamental conclusions. In the frameworks of 
2D calculations, a quasi periodic consequence of "on/off"' phases 
is expected. For the computational B -field strength the typical pe- 
riod of these oscillations is a fraction of a second. In fact for a large 
range of the wind magnetic field value, this period depends only on 
the accretion rate and BH mass (see for detail Sect.|3]). 




Figure 4. The flow configuration corresponding to the "magnetically ar- 
rested" regime. The solid lines show the magnetic field lines and the arrows 
indicate the flow velocity (both the direction and the absolute value). The 
color corresponds to log p. It can be seen that the accretion is stopped at 
this moment and jet ejection is significantly less efficient as compared to 
the configuration shown in Fig. [2] The initial magnetic field in the stellar 
wind was assumed to be 10^ G, and the accretion rate was selected to be 
Mbh = 8 X 10~^^Moyr~^. The shown snapshot of the flow corresponds to 
the time moment f = 0.159 s after the accretion was initiated (this moment 
corresponds to the "paused jet" as seen in Fig.[3}. 

3 DISCUSSION 

The obtained numerical solution may have some interesting im- 
plications for the phenomenological study of yuQs. Namely, it was 
shown that an efl'ective jet launch is possible in binary systems in 
the case of the direct wind accretion, i.e. without the formation of 
an accretion disk. Since the present phenomenological description 
of the processes occurring in yuQs relays on the observational ap- 
pearance of the accretion disk, which lacks in the case of direct 
wind accretion /iQs should fall out from this framework. Although 
the in-falling matter gets strongly heated at the termination shock 
of the in-falling matter, given the very short free-fall and dynami- 
cal time- scales compared to the cooling time, no detectable thermal 
emission is produced in the framework of the scenario. Thus, the 
only way to obtain observable features of the suggested scenario 
is to study in detail the properties of the jet, and in particular its 
modulation properties. 



3.1 Jet internal time-scales 

As it was discussed above the reason behind the strong modulation 
("on/ofl"' states) of the jet power is the accumulation of the mag- 
netic field flux at the BH horizon. When the magnetic field pressure 
exceeds the ram pressure of the accretion flow, accretion is stopped 
and the jet switches ofl". It is possible to estimate the variability time 
from the comparison of the two key forces: magnetic pressure 

p - - n\ 
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Figure 1. The square root of the jet luminosity as a function of the magnetic flux, ^, accumulated at the BH horizon. Three diff'erent regimes can be seen in 
the figure: (i) the initial accumulation of the magnetic flux (marked with "1"); (ii) transition to the BZ regime (marked with "2"); and (iii) the conventional BZ 
regime, when a linear dependence L?^^ oc Y is expected (marked with "3"). The initial magnetic field in the stellar wind was assumed to be 10^ G, and two 
accretion rates were considered Mbh = 8 x 10~^^Moyr~^ (left panel); and Mbh = 1.5 x 10~^^Moyr~^ (right panel). 



and gravitational hydrostatic pressure 



_ T GMbhM 



(8) 



Here A 1 is a dimensionless parameter which accounts for the 
flow geometry; r is the time since the moment when accretion is 
stopped, and G is the gravitational constant. Combining these two 
relations one obtains the time required to restart accretion: 



VJ/2 



2;r2GMBH M 



(9) 



If there is no escape of the magnetic flux, i.e. the value of ^ is 
increasing monotonically with time t, the interval between jet acti- 
vations T becomes longer with time, and asymptotically, given the 
quadratic dependence on ^ in Eq.(|9|, the accretion process stops 
completely. Thus, the key condition for the quasi-periodic jet is the 
escape of the magnetic flux from a magnetically arrested region. 
We note that the escape of the magnetic field can be studied numer- 
ically only in the a 3D approach, i.e. a self-consistent description 
of this process remains out of the scope of this paper, and should 
be studied elsewhere. In what follows we adopt the critical value 
of ^cr. when the flux escapes, as a free parameter. This parameter 
is in fact directly connected to the jet power through Eq.([6|, thus 
the escape of the magnetic field defines not only the duration of the 
interval when the jet is paused, but also its available power during 
the jet ejection. 

In the case of non- steady jet activity, the averaged jet lumi- 
nosity is the key quantity characterizing the possible observational 
appearance of such jets. Obviously, the averaged power may de- 
pend on many dififerent aspects of the involved processes, e.g. on 
the magnetic flux escape rate. In the simplest case, the critical mag- 
netic flux value is smaller than the value required to pause the ac- 
cretion: < i-e. the magnetic flux starts to escape before the 
magnetic field can stop the accretion. Then, the jet is expected to 
be persistent with a power of 



J^BZ ~ ^cr,20 



Mbh 
3Mo 



erg s 



(10) 



where ^cr,20 = ^cr/lO^^Gcm^. Given the assumed condition < 
an upper limit for the jet luminosity can be obtained. Our nu- 
merical simulations show that the jet power Lj = CMc^, where 
C ~ 0.05, corresponds to the magnetic flux Here we note that, 
contrary to the widely accepted disk accretion paradigm, in the con- 
sidered case the accretion rate cannot be estimated from the thermal 
radiation produced by the accreted matter. 

In the case if the > one should expect a more compli- 
cated temporary patterns in the system. Namely, the quasi-periodic 
behavior of the jet will not be vanished by the escape of the mag- 
netic field, and the jet active phases will be separated by period of 
the jet silence of duration 



VJ/2 



(11) 



Although the interval of the jet active phase depends on the process 
of the magnetic flux escape rate, the averaged jet luminosity can be 
estimated assuming that the escape process is rather slow as com- 
pared to the variability time. In this case, the duration of the active 
phase is roughly determined by time interval during which the ac- 
cretion flow can keep the magnetic flux of near the BH horizon. 
Accounting for the phases when the accretion stops, the matter ram 
pressure can be estimated as 



Pram ~ 



47tR] 



,5/2 



(12) 



where fac is the duration of the jet active phase, and Rq is the gravi- 
tational radius of the BH. Using Eq.(|7]), one can obtain the follow- 
ing simple relation 



^cr ~^ ^ac 



(13) 



^ We note that the value of the parameter C depends on the BH rotation 
parameter a. Since the critical value of the magnetic flux ^cr is a weak 
function of a (see e.g. |Tchekhovskoy et al.|20lT} , the parameter C should 
be roughly proportional to the function f{a) from Eq.(|6}. Thus, for a fast 
rotating BH, this parameter can have a value up to 10 times larger than the 
one obtained in our simulations. 
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Figure 2. The magnetohydrodynamical properties of the flow: the solid lines show the magnetic field lines and arrows indicate the flow velocity (both the 
direction and the absolute value). The top panels show the distribution of the mass density (logp is shown by color) for the whole computation regions (top 
right panel), and for the zoomed central part (top left panel). The bottom panels indicate the ratio of the toroidal to poloidal magnetic field components 
(log Bfp/ Bp is shown by color in the bottom left panel); and the ratio of the gas to magnetic field pressures (logPg/Pm is shown by color in the bottom right 
panel). The initial magnetic field in the stellar wind was assumed to be 10-^ G, and the accretion rate was selected to be Mbh = 8 x 10~^^Moyr~^ . The shown 
snapshot of the flow corresponds to the time moment r = 0.145 s after the accretion was initiated. 



Since the BZ luminosity is proportional to the jet averaged lu- 
minosity appears to be independent on the ^cr- 



<Lj 



'^cr "I" ^ac 



XLBz(^cr)=Z^Bz(^0)^CMc2 



(14) 



3.2 LS 5039as a jiQ 

The universal averaged jet power expected in the suggested sce- 
nario allows to make meaningful estimates for the specific case of 
LS 5039 without 3D modeling. Indeed, the jet power depends on 
the accretion only and for the expected rate of M = 8 x 1 ^ Moyr~^ 
one obtains 



Thus, for a broad range of parameters one can expect a universal 
dependence of the jet power on the accretion rate. 



Lj,Ls ^3x10^^ ergs- 



CIS) 
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Figure 3. The jet luminosity as a function of time, since the moment when the accretion was started. The initial magnetic field in the stellar wind was 
assumed to be 10^ G, and two accretion rates were considered Mbh = 8 x 10~^^Moyr~^ (left panel); and Mbh = 1.5 x 10~^^Moyr~^ (right panel). 



per each jets. This estimate is similar to flux level detected in GeV 
energy band with FermifLAI telescope, thus to explain total non- 
thermal luminosity in the framework of the //Qscenario one needs 
to invoke Doppler boosting of the emission. The Doppler boosting 
factor depends on the jet bulk Lorentz factor and the orientation of 
the jet in respect to the observer. The jet is orientated along the BH 
rotation axis, thus most likely it is directed perpendicularly to the 
orbital plane. The most feasible orbital inclination in LS 5039 is 
/ ~ 20° - 25° ( [Casares et al.|2005|[S^ty et aL|201l] ), thus one can 
obtain an upper limit on the boosting factor: S = l/[r(l-/3 cos /)] ^ 
2r/(l + r^f). If r > l// the emission emitted toward Earth will 
be de-boosted, the optimal boosting can be achieved if F ^ 1 //. 
Thus, the largest Doppler boosting factor can be estimated sls S < 
3, which should lead to the apparent jet luminosity at the level of 
< lO^^ergs"^ For a relatively high, > 10% conversion rate of the 
jet luminosity to non-thermal emission, this energetics can explain 
the observed broad-band emission, although a dedicated modeling 
is required to shed the light on the possibility of the production of 
the detected radiation in the framework of the fiQ scenario. 



4 SUMMARY AND FINAL REMARKS 

In this paper we study the wind accretion onto a rotating BH in 
close binary systems harboring a young massive star. It is shown 
that the angular momentum of the accreted stellar wind material is 
not sufficient for the formation of the accretion disk. On the other 
hand, the direct wind accretion can create conditions for the activa- 
tion of the BZ process, and a powerful jet can be launched without 
any signatures of accretion, unlike it is usually seen in bright ther- 
mal X-ray binaries. The jet power is expected to be modulated on 
a sub-second scale, which can lead to an efficient jet loading with 
stellar wind material (see Derishev et al. 2011, in preparation). Al- 
though the flux escape mechanism should be properly treated with 
regard to the formation of persistent/quasi-periodic jets, a universal 
estimate on the jet averaged power is obtained. It is expected that 
approximately ~ 10% of the accreted rest energy can be chanelled 
into the jet. In the specific case of LS 5039, the obtained jet lumi- 
nosity can be responsible for the observed non-thermal radiation 
if one invokes Doppler boosting, which can enhanced the apparent 
flux from the system and to be additional factor of the modula- 



tion during the orbital motion. The future study of the non-thermal 
production in the framework of the suggested scenario, would re- 
quire a self-consistent modelings of the wind accretion to the BH 
( [Owocki et al.|2011|), propagation of jet thro ugh binary system en- 
vironment ( |Perucho & Bosc h-Ramo n|2008[[Perucho et al.|2010| ), 
and non-thermal processes occurring in the jet. 
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